arXiv:1507.04872v2 [cond-mat.quant-gas] 18 Dec 2015 


Hybrid synchronization in coupled ultracold atomic gases 


Haibo Roberta Zambrini,^ Artur Polls,Joan Martorell,^ and Bruno Julia-Diaz^’® 

^Departament d’Estructura i Constituents de la Materia, 

Universitat de Barcelona, 08028 Barcelona, Spain 
^College of Science, Xi’an University of Posts and Telecommunications, 710121 Xi’an, China 
^Institut de Ciencies del Cosmos, Universitat de Barcelona, 
lEEC-UB, Marti i Franques 1, E-08028 Barcelona, Spain 
’^Instituto de Fisica Interdisciplinar y Sistemas Complejos IFISC (CSIC-UIB), 

Campus Universitat Hies Balears, E-07122 Palma de Mallorca, Spain 
^ICFO-Institut de Ciencies Fotoniques, Parc Mediterrani de la Tecnologia, 08860 Barcelona, Spain 

(Dated: December 21, 2015) 

We study the time evolution of two coupled many-body quantum systems one of which is assumed 
to be Bose condensed. Specifically, we consider two ultracold atomic clouds each populating two 
localized single-particle states, i.e. a two-component bosonic Josephson junction. The cold atom 
cloud can retain its coherence when coupled to the condensate and displays synchronization with 
the latter, differing from usual entrainment. We term this effect among the ultracold and the 
condensed clouds as hybrid synchronization. The onset of synchronization, which we observe in the 
evolution of average properties of both gases when increasing their coupling, is found to be related 
to the many-body properties of the quantum gas, e.g., condensed fraction, quantum fluctuations 
of the particle number differences. We discuss the effects of different initial preparations, and the 
influence of unequal particle numbers for the two clouds, and we explore the dependence on the 
initial quantum state, e.g. coherent state, squeezed state, and Fock state, finding essentially the 
same phenomenology in all cases. 

PACS numbers: 


I. INTRODUCTION 

Synchronization has been described in physics, chem¬ 
istry, biology and social behavior m- It has been ex¬ 
tensively studied in classical non-linear dynamical sys¬ 
tems |3], and chaotic ones [5]. The same phenomena 
have been explored recently in quantum systems, e.g. 
opto-mechanical devices [7], damped harmonic oscilla¬ 
tors [6l [TT], driven [8] and purely dissipative spins [9], 
and non-linear optical cavities uni. Synchronization can 
refer to the mutual effect between detuned but otherwise 
equivalent components adjusting their rhythms (sponta¬ 
neous synchronization) as, for instance in Refs. [H IHl E] ■ 
Otherwise, a slave system can be driven to follow the dy¬ 
namics of an external source leading to entrainment or 
driven synchronization, as for instance in Refs. [7l |8] . In 
quantum many-body physics connections between quan¬ 
tum entanglement and mutual synchronization have been 
discussed in continuous variable systems Enmna. 

Ultracold atomic gases are particularly relevant quan¬ 
tum many-body systems. Since the first experimental 
production of Bose-Einstein condensates (BEC’s), they 
have evolved from being a theoretical curiosity to ver¬ 
satile systems potentially useful in a large number of 
fields [13) . Identifying the onset of synchronization in 
these systems and proposing ways in which such phe¬ 
nomena can be characterized both experimentally and 
theoretically is a significant step forward in our under¬ 
standing of the dynamical evolution of coupled quantum 
many-body systems. 

Among the most promising applications are those that 
stem from the macroscopic sizes of the condensates. 


BECs are fantastic candidates for high accuracy inter¬ 
ferometric devices [HI [15]. These devices rely on the 
high degree of coherence maintained by BECs. 

In recent years, experiments with bi-modal ultracold 
atomic gases have managed to produce entangled en¬ 
sembles in which the interferometric capabilities can be 
largely enhanced [T6l - f20j . These improved interferomet¬ 
ric properties are directly related to the pseudo-spin 
squeezing which can be produced using several tech¬ 
niques [2TH2^ . One of the main sources of decoherence in 
BECs are atom-atom interactions which induce dephas¬ 
ing of the different Fock components [H]- Several pos¬ 
sibilities, notably the generation of squeezed states [23], 
have been proposed to increase the coherence times and 
improve the interferometric capabilities, e.g. of the re¬ 
cent Mach-Zehnder proposal [25) . 

In this paper we describe how decoherence effects due 
to the atom-atom interaction can be largely suppressed 
if a quantum many-body system is coupled to a Bose- 
Einstein condensate. To be more specific, we consider 
two bosonic Josephson junctions, a and b. Subsystem 
b is taken to be a BEC at all times during the evolu¬ 
tion [26) . while subsystem a is a standard bosonic Joseph¬ 
son junction, i.e. it may fragment during the evolu¬ 
tion [2HI EH EH] • The coupling between a and b is pro¬ 
vided by the contact inter-species interaction. Directly 
related to this improvement in the coherence times of 
subsystem a, is the onset of a hybrid synchronization be¬ 
tween a and b. 

This synchronization is called hybrid because the two 
coupled atomic samples are in different regimes, one be¬ 
ing a BEC described within a classical approximation and 
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the other being a cloud of cold atoms described with a 
fully quantum formalism. We notice that in the literature 
the term hybrid synchronization has also been used in 
other contexts to describe, for instance, synchronization 
between (classical) chaotic systems persisting despite the 
difference in some variables [29] or to the co-existence of 
different kinds of synchronization in composed (classical) 
systems like cascade-coupled lasers Ref. m- 

The paper is organized as follows: In Sect. |ll| we in¬ 
troduce the model system, a two species two-site Bose- 
Hubbard model. Assuming that one of the species is 
condensed, we construct our mixed quantum-classical de¬ 
scription. In Sec. the coupled dynamics is studied, 
particularly the onset of hybrid synchronization and its 
relation to the condensed fraction of subsystem a. Fi¬ 
nally, in Sec. IV we present our conclusions. 


aa, bb and ab contact interactions. The coupling between 
the two gases is solely due to the term proportional to 
Uab, which is an on site interaction between the atoms 
of the two species. As customary, we introduce the def¬ 
inition of the population imbalance of each species as, 
Za = {na,L - na,R)/Na and Zb = {hb,L - nb,R)/Nb. 

The quantum [IT] and classical (fully condensed) dy¬ 
namics [T^HTS] , of this model have been previously stud¬ 
ied comparing different dimensional reductions of the 
three-dimensional equations |47j . Also, measure syn¬ 
chronization [35] has been studied in the classical [35] 
and full quantum case [50] . 


A. Hybrid, quantum and classical, description of 
the system 


II. MODEL 

Let us consider two different kinds of bosons, a and 
b, populating two single-particle states, L and R. The 
main results described in this paper are essentially inde¬ 
pendent of the nature of the single-particle states. A way 
to produce them is by trapping ultracold atoms in a spa¬ 
tial double-well potential [IHISI]. A different possibility 
would be to populate two hyperfine states of the atom 
and couple them linearly as in Ref. [32] . A solid-state al¬ 
ternative could be provided by extending the experiments 
in Ref. [33] to two different kinds of exciton-polaritons. 

The interaction between the bosons is assumed to be 
described by a contact term with strength proportional to 
the s—wave scattering length. For simplicity the bosons 
are taken to be two different hyperfine states of the same 
atom, so that the masses of both species a and 6, is the 
same. With these simplifications, the many-body, Bose- 
Hubbard, Hamiltonian for the system reads, 

n = 'Ha + nb + 'Hab ( 1 ) 


where 


H 


a 


Hb 


idab 


-Ja[d\dR -I- OLaJj) 

^{na,L{na,L “ 1 ) + na,R{ha,R - 1 )), 

-MblhR + hV^) 

{nb,L{nb^L — 1 ) + nb^R{hb^R — 1 )), 

^ («a.L - na,R){hb,L “ fjfc.fl) ■ (2) 


^:,{R)(^L{R)) are creation (annihilation) 
operators for the L or R modes of a and b. The Hamil¬ 
tonian includes tunneling terms, proportional to Ja(b), 
which the in absence of any interaction induce peri¬ 
odic Rabi oscillations of the populations between the 
states. The contact interaction translates into terms with 
strength proportional to Ua, Ub, and Uab which gauge the 


Our main interest is to study the combined evolution 
of a and b. In particular, one of our aims is to discern 
whether the coupling to a condensed system, 5, will en¬ 
hance the degree of condensation of the ultracold gas, a. 
To answer this question, b is assumed to be condensed 
at all times. Physically, this situation could be attained 
if the b component has a large enough number of atoms. 
We thus neglect the quantum fluctuations of the b cloud. 
In practice this is done by replacing the operators by 
c—numbers, br = , bR = We de¬ 

fine the phase difference as = (pb^R — 4>b,L- 

The Hamiltonian is thus written as, H = Ra + Hb + 
Hab, with 

Hb — ‘^Jby/blb^Rnb^R COS (pb 

+ {nb,L{nb,L — ^) + nb^Rim^R — 1)) , 

Hab = ^ina,L-na,R)iNbZb). (3) 

where Hb is the c—number version of Hb- To study the 
time evolution we solve the following coupled set of equa¬ 
tions, 

ldt\'i!a) = H{t)\^a), (4) 

Zb = -2Jb^l - zl sincpb, 

pb — ‘^Jb-^bZbR 2Jb . == COS (pb T 2Jb-^ab{Za') (t') ■ 

V 1 “ h 

where we have introduced the dimensionless ratios, A^ = 
NaUa/i2Ja), Af, = NbUb/{2Jb), and Kb = NaUab/{2Jb). 

The time dependence of H stems from the time depen¬ 
dence of Zb and (pb- This set of equations is obtained con¬ 
sidering the b subsystem as a classical parametric driving 
for the a subsystem and including the feedback effect of a 
on b. Neglecting quantum fluctuations in one field (then 
classical), still retaining feedback effects due to the in¬ 
teraction with the other (fully quantum), is reminiscent 
of the ‘time-dependent parametric approximation’ [33], 
used to describe large quantum fluctuations in a convec- 
tively unstable signal of an optical parametric oscillator. 
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FIG. 1: (Color online) Evolution of the distribution of Fock coefficients, |cfep in Eq. ([^ for the quantum gas a as a function 
of time. We show three different values of the coupling A^f, = 0 (a), A^f, = 0.1 (b), and A^f, = 1.0 (c). The insets depict the 
distribution \ck\^ (linked by a line) for three specific times t = tuabi, Stuabi, and ISfRabi and the three different couplings. 
In all cases Na = Nf, — 30, Ja = Jt = 1 and A^ = A;, = 1. The initial states are described in the text, (za) = Zb = 0.4 
and (jia = (jib = 0. The value of \ch\^ (color) is plotted as a function of the imbalance of each Fock state, Za,k = {2k/Na — 1), 
fc = 0,..., Na- 


when the pump field is approximated by a classical field 
while the signal is quantum |35j . Furthermore, here we 
are neglecting the effect of the quantum fluctuations of 
Za in evaluating the dynamical evolution of (j)b, as we 
approximate Za by its expectation value. 

The coupled system of Eqs. Q is solved in the follow¬ 
ing way. We use a fourth-order Runge-Kutta routine to 
integrate the differential equation for zj, and coupled 
to a unitary truncation of the Schrodinger equation for 
Idfa) of the form. 


|«'a(t+At)) 






l»4<)) 

(5) 


with At « 0.0002 tRabi- 

The solutions of Eq. Q are numerically found to con¬ 
serve the average energy, (T'a|iT(t)|5'a) ~ (Ta|.ff(0)|T'a) 
in all the calculations reported in this paper. It is worth 
emphasizing that the coupling between a and b cannot 
simply be regarded as a driving term for a. Conservation 


of the total average energy implies on average a transfer 
of energy between the a and b subsystems. The dynamics 
is thus radically different from the case of a driven sin¬ 
gle component Josephson junction [Si [37], which for in¬ 
stance would occur if Zb{t) in our description was replaced 
by a periodic function. In such case, energy would not 
be conserved and subsystem a would gradually increase 
its energy. In our formulation the exchange of energy al¬ 
lows for the mutual synchronization of both species and 
differs from entrainment as will be described in the fol¬ 
lowing section. 

Using the Fock basis of the Na + 1 dimensional space, 
\na,L,na,R) = {|lVa,0),..., |0,iVa)}, the most general a 
state is written as, 


Na 

I'ita) ='^Ck \k,Na - k). (6) 

k=0 

All many-body properties of the state are computed from 
the CfeS, e.g. the average population imbalance of the a 



















4 





0 5 10 15 20 25 30 35 40 45 

^'^Rabi 


FIG. 2: (Color online) Evolution of two relevant many-body 
properties for different values of the interspecies coupling, 
Aai,. Panels (a) and (b) depict the dispersion of the popula¬ 
tion imbalance of the a cloud, and the condensed fraction 
of the a system, ni, respectively, for the values considered in 
Fig. Aai, = 0, Aab = 0.1, and Aoi, = 1.0. Panel (c) scru¬ 
tinizes the evolution of the condensed fraction for a broader 
range of values of 0 < Aai, < 1. All other parameters are the 
same as in Fig. Note that the initial condensed fraction, 
ni( 0 ) = 1 , so that the plotted values are rii. 


cloud reads. 


No. 

i^a) — ^ ^ \^k\ ^a,k 
fc=0 


(7) 


with Za,k = {2k/Na — 1). The degree of condensation 
of the cloud is given by its condensed fraction, i.e. the 
largest eigenvalue, ni, of the single-particle density ma¬ 
trix p“(t), pfj{t) = ;^(4'a(t)|aJaj|4'a(i)): bJ = L,R. 
The condensed fraction is also referred to as single¬ 
particle coherence. With this normalization, a fully con¬ 
densed cloud of a would correspond to ni = 1 and n 2 = 0. 
In the following sections we discuss the time evolution of 
ni for different initial conditions and couplings. 



FIG. 3: (Golor online) Time averaged condensed fraction ni = 
fg ni(t)dt with T = IStRabi as a function of Aai, and Aa = 
A;,. The initial condition is again fully condensed for the a 
component with ( 2 a) (0) = 25 (0) = 0.4. All other parameters 
are as in Fig. The dashed red lines correspond to 1 -|- Aa — 
Aai, = 0 and 1 -|- Aa -|- Aai, = 0, see text for details. Notice 
that the scale in the color bar is nonlinear, and it spans all 
the values of the data. 


III. COUPLED QUANTUM AND CONDENSED 
DYNAMICS 

In this section we will describe how the appearance 
of synchronization in the combined evolution is found 
to be directly related to a coherent, in the sense of not 
dephased, evolution for the a subsystem. 


A. From dephasing to coherent evolution 

Let us first exemplify our discussion with one spe¬ 
cific configuration. We will choose as initial state a 
condensed quantum state for a, i.e. all a atoms popu¬ 
late the same single particle state l/-\/2(cos(61a/2)a^ -|- 
g*0a sin(61a/2)a^)|vac). The many-body state reads, 

/A 

= E [cOs(da/2)]'=[sin(0a/2)]^-'= 

k^O ^ 

X . (8) 

In particular we will take cos(da) = (2a}(i = 0) = 0.4 
and (j)a = 0. To emphasize the effect of the coupling 
term we choose similar conditions for b, Zb = 0.4, (j)b = 0, 
Na = Nb = 30, J = Ja = Jb = i and = Af, = 1. 
The Rabi time is tRabi = tt/J, and the Rabi frequency is 
a^Rabi — 27r/tRabi = 2T 

In absence of coupling between a and b, Aab = 0, 
the quantum system a with a non-zero initial popula¬ 
tion imbalance evolves with time in a well studied fash- 
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FIG. 4: (Color online) Evolntion of {za){t) (dashed-red), and 
Zi,{t) (solid-black) for two different values of the coupling (a) 
Aab = 0 and (b) Aab = 0.1. All other parameters are as in 
Fig.[T] 

ion [351 |3H1 ISn] ■ Due to the atom-atom interactions which 
dephase the different Fock components, the initial distri¬ 
bution of Cfc evolves in time deforming its initial shape 
(see Fig. I^a)). For the first oscillations, up to t ~ StRabi 
the wave packet remains mostly unchanged, which in turn 
is also reflected in the fact that the a component remains 
essentially condensed rii > 0.98 (see Fig. ib)). For 
larger times, t > lOtRabi, the original shape is lost, the 
a component is no-longer in a coherent quantum state, 
and thus the condensed fraction drops below 0.95. Inter¬ 
action among atoms a is thus seen to decrease the degree 
of condensation of the subsystem fairly early. Directly 
related is the increase in the uncertainty on the particle 
number difference, erf = (zf) — (zq)^, shown in Fig. [^a), 
also appreciable in Fig. [^a). 

When coupling the quantum system a to the condensed 
one, 6, a distinctive dynamics is found. The dephas¬ 
ing due to the atom-atom interaction disappears and the 
quantum system remains condensed for longer times. As 
seen in Fig. Bb) for Aab = 0.1 the distribution of |cfcp 
remains closer to a displaced binomial one, which again 
reflects in a much larger condensed fraction (see Fig. 
and smaller insets in Fig. [^. Already with this fairly 
small value of Aab we find a substantial increase in the 
condensed fraction, which is now at all times larger than 
0.97. Further increasing Aab, the effects are enhanced: 
The quantum system remains close to condensed for long 
times, see Fig. [^c), and the distribution of |cfcp thus 
evolves, keeping its origmal shape. The latter is also 
reflected in erf, see Fig. |^a), which is found to remain 
almost constant in the coupled case. Effectively, the cou¬ 
pling to the condensed gas removes the dephasing effects 
due to the atom-atom interaction, obtaining an almost 
interaction-free evolution of the quantum system. 

This phenomenon is persistent in a broad range of pa¬ 
rameters. For Aab 0.01 and maintaining similar initial 
conditions, the quantum evolution is essentially coherent, 
as shown in Fig. (b) and [^c). 


Varying the values of Aa = Ab and Aab one finds the 
following picture. One can easily prove that within our 
formalism for Aa = 0, irrespective of the value of Aab 
the condensed fraction of a remains 1. This behavior 
survives for small values of Aa < 0.5, for which a high 
degree of condensation for a is also found (see Fig. |^. 
For Aa > 0.5, increasing Aab tbe condensed fraction of 
the system is found to remain essentially constant in time 
(see Fig.|^. Thus, the coupling between a and b increases 
the coherence of the a system. However, for larger val¬ 
ues of Aab, in particular for Aab > Aa -I- 1 we observe a 
decrease of the condensation of a (see Fig. [^. This can 
be understood from the linear stability analysis of the 
classical equations around Za,Zb I. In this case, one 
of the two natural modes m ^2 — WRabiVl + Aa — Aab 
becomes unstable if Aab > Aa -I- 1, which induces deco¬ 
herence of the a system. 

A similar picture is obtained in the attractive inter¬ 
species interaction case, Aab < 0 (see Fig. |^. In particu¬ 
lar, for values of Aa < 4 we see that increasing the value 
of I Aab I the condensation of the a cloud is increased. Fur¬ 
ther increasing |Aab|, as observed for Aa ~ 1.5, the degree 
of condensation decreases. The boundary of the classical 
stability region obtained for small values of the imbal¬ 
ance is found by imposing wi = WRabiVl + Aa -I- Aab to 
be real, Aab > -Aa - I. 


B. Hybrid synchronization 

We have described how the coupling between the sub¬ 
systems prevents the a subsystem from fragmenting dur¬ 
ing the time evolution for certain coupling values. Now 
we show how this effect is directly connected to the ap¬ 
pearance of synchronization between properties of both 
subsystems. This synchronization, to which we refer as 
hybrid, stemming from the hybrid nature of our coupled 
system, manifests itself in average properties. 

In the non-coupled case. Fig. (a), the population 
imbalance of the condensed subsystem, Zb, is fully peri¬ 
odic [35]. The frequency seen in the figure, ui ~ l.SSwRabi 
is close to the one obtained linearizing around the Zb = 0 
fixed point, uj = WRabi vT+^Ab ~ WRabi V^. The popula¬ 
tion imbalance of subsystem a, {Za), features character¬ 
istic collapses and revivals, which are also present in the 
condensed fraction and population imbalance dispersion 
shown in Fig. [38]. 

Coupling both subsystems. Fig. ^(b), both signals are 
found to be much more correlated. To quantify this we 
compare in Fig. [^the frequency spectra of these signals 
for different values of Aab- In the uncoupled case the 
quantum signal is found to have several peaks around 
the same frequency lo ~ 1.38a;Rabi. The different equis- 
paced peaks reflect the long evolvent seen in Fig. (upper 
panel). They arise from the atom-atom interaction which 
makes the spectrum of the many-body Hamiltonian de¬ 
part from the equispaced/harmonic case in the uncoupled 
case producing quantum revivals [31110]. 


























































































































































6 





FIG. 5: (Color online) Absolute value of the frequency spec¬ 
trum z(uj) (Fourier Transform of the average population im¬ 
balances), for different values of the coupling Aab- The fre¬ 
quency spectrum of the quantum result, {za}{ijj) is given 
in the left panels. Zb(oj) is depicted in the right panels. 
The frequency spectra are obtained from time series up to 
Tlviax “ 585 tuabi. 


^^Rabi 


FIG. 6: (Color online) Time correlation coefficient 
defined in Eq. versus the coupling strength Aai,. Panel (a) 
corresponds to two specific values of Aai, = 0, Aab ~ 0.1, while 
(b) contains the time evolution of ^ broader set 

of Aab- At = 14tiiabi- All Other parameters are as in Fig. 


As Aab is increased, the spread of the peaks in the 
quantum case is reduced. For Aab = 0.1 the Fourier de¬ 
compositions of both signals are very similar, showing 
a large peak at wtRabi — 1.405. For this value of Aab, 
a is mostly condensed and the classical description of 
the full system should approximately hold [HJ 03] . In¬ 
deed, the found frequency is reproduced by the classical 
equations [43]. For this particular case of similar initial 
conditions of a and b, the classical description of the bi¬ 
nary mixture predicts, linearizing around z <C 1, just one 
frequency, = WRabiVl + -b Aab = V^WRabi [43]. 
The deviation observed is due to the departure from 
z ^ 1 of our initial conditions (2:a)(0) = Zb{0) = 0.4. 

The synchronization phenomenon, as seen in the 
Fourier analysis of Fig. goes from several different fre¬ 
quencies for both subsystems in the uncoupled case, to 
a single major frequency in the coupled case. Thus the 
two subsystems get frequency-locked as the interaction 
is increased. It is also worth noting again that in our 
description subsystem a can also remain condensed, but 
b is not allowed to fragment. The resulting scenario is 
that the coupling induces condensation in subsystem a. 

In order to have a quantitative characterization of the 
synchronization, we calculate the time correlation coef¬ 
ficient C, which can be used to judge whether two time 
series are synchronized um- For two time signals, /(t) 
and g(t), it is defined as. 


Sf{t)Sg{t) 
5p{t) Sg^{t) 


where the bar stands for a time average f{t) = 
S dt' f{t') with time window At and Sf{t) = 

f{t) — f{t). We choose a time window At = 14tRabi. 
For in phase (anti-phase) synchronization C ~ 1 (—1), 
while it equals zero for fully non-synchronized cases. 

In Fig. 1^ we show C as a function of time for dif¬ 
ferent values of the coupling strength Aab. The figure 
has a similar structure as that found when computing 
the condensed fraction of the same subsystem in Fig. 
The subsystem is seen to remain condensed for the same 
values of Aab for which the two signals are synchronized. 

For Aab < 0, we have shown that the results were sim¬ 
ilar to the Aab > 0 case, that is the condensation for the 
a cloud is enhanced with enough coupling (see Fig. |^. 
As occurred in the Aab > 0 case, the increase in conden¬ 
sation for Aab < 0 is also accompanied by an increase in 
the time correlation function between the average popu¬ 
lations of both species, e.g. for Aab = —TO with the same 
parameters and initial conditions as in Fig. [^ C ~ 1 is 
reached. 


1. Effect of unequal populations 

The hybrid system we are considering, with b remain¬ 
ing fully coherent during the time evolution, is justified 
if the number of b atoms is large enough. Up to now we 
have discussed the case in which Na = Nb, in order to 
isolate the effect of the inter-species coupling. 

Now, we consider a system with unequal population 
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FIG. 7: (Color online) Evolution of single particle coherence 
ni for different values of Aab- With unequal population ^ = 
10. In all cases we fixed At, = 1, A;, = lOAa, and Nb = 300, 
Na = 30. The initial state and other parameters are as Fig.[2 
(a) corresponds to the two values of Aab- (b) scrutinizes the 
same function for a more detailed range of values of 0 < Aab < 
1 . 


Nb/Na = 10. With Nb = 300, = 30, and A;, = lOA,,, 

Aa = 1, so that we have Ua = Ub- We keep the initial 
state and other parameters the same as in Fig. The 
results show a picture very similar to that of the case with 
equal populations. In Fig. we show the evolution of 
single-particle coherence, ui, for different values of A^f,. 
We observe an increase of the condensed fraction as Aab 
is increased. This is similar to what has been shown in 
Fig.H, However, to make a more quantitative comparison 
of Fig. [^c) and Fig.[^(b), we will notice that for the case 
with unequal populations, in order to reach the same level 
of single particle coherence as for the case with Na = Nb: 
a relatively larger Aab is needed now. 

Fig.lH shows the hybrid synchronization for this case 
[Nb ^ Na). In the non-coupled case, Fig. (a), the 
population imbalance of the condensed subsystem, Zb, is 
fully periodic. Compared with the equal population case, 
we notice that the frequencies of the two signals are very 
different. However, by coupling both subsystems with 
Aab = 0.1, as seen in Fig. (b), the signals are found 
to be more correlated. Further increasing Aab, the two 
signals show synchronous dynamics due to this coupling 
effect [Fig. (c)]. 

To quantify the hybrid synchronization, in Fig. we 
show the evolution of C(^za),zb> for different values of Aab- 
With Aab = 0, Fig. [^a), the time correlation function is 
close to zero as expected. With Aab = 1, Fig. [^a), the 
time correlation function is close to I as hybrid synchro¬ 
nization occurs. By comparing Fig. |^b) with Fig. j^b). 


FIG. 8: (Golor online) Evolution of {za){t) (dashed-red), 
and Zbit) (solid-black) with Ab = lOAa, A^ = 1, and unequal 
populations Na = 30, Nb = 300. Three different values of the 
coupling are chosen, with (a) Aab = 0, (b) Aab = 0.1, and (c) 
Aab = 0.6. All other parameters are as in Fig. 


one sees that in order to reach hybrid synchronization 
in the case with unequal populations, a larger coupling 
strength Aab is needed. This is a general feature of syn¬ 
chronization arising when the coupling between systems 
is large enough to overcome their detuning. Indeed here 
the detuning between the two clouds increases with the 
difference between the populations and needs to be com¬ 
pensated by a larger reciprocal contact interaction. 


2. Effects of different initial preparation 


Up to now we have considered the same initial popu¬ 
lation imbalance for a and b. Now we consider a more 
general case, in which the initial average population im¬ 
balance is not the same for both subsystems. In partic¬ 
ular we will consider (za) = 0.4 and Zb = 0.2, while all 
other conditions are as in Fig. The main finding de¬ 
scribed earlier is again found: Coupling a and b increases 
the coherence time of subsystem a. As described above, 
the fragmentation of subsystem a which takes place dur¬ 


ing the uncoupled evolution, see Fig. decreases as 
the coupling between a and b is increased. For values of 
Aab — Aa = Ab, a remains almost fully condensed during 
the evolution. 

Let us analyze the Fourier decomposition of the evolu¬ 
tion of the population imbalance. Following our previous 
discussion in the case of decoupled clouds, we find that 
Zb has only one frequency, which now is closer to the 
one expected in the linear regime, w ~ WRabiVl + Af, = 
■\/2wRabi- The a subsystem, as before has a number of 








































FIG. 9: (Color online) Evolution of for different val¬ 

ues of Aab- With unequal population ^ = 10. In all cases we 
fixed Aa = 1, A;, = lOAa, and Nb = 300, Na = 30. The initial 
state and other parameters are as Fig. (a) corresponds to 
the two values of Aab- (b) scrutinizes the same function for a 
more detailed range of values of 0 < Aab < 1. 


FIG. 10: (Color online) Evolution of the condensed fraction 
of subsystem a, ni, as a function of time for different values of 
the coupling Aab- The upper panel corresponds to the three 
values, Aaf, = 0, 0.1, and 1.0. The lower panel scrutinizes the 
same function for a more detailed range of values of 0 < Aab < 
1- Here we consider different initial conditions, (za)(0) = 0.4, 
and 2 :i,( 0 ) = 0.2. 


peaks, whose spread is related to the deviation from the 
equispaced spectrum (which in this case is smaller). As 
we couple the two subsystems, the spread in the a sub¬ 
system disappears and two prominent frequencies appear 
both for a and b- Already for A^t = 0.45, we have the 
same Fourier structure in both cases, signaling the ap¬ 
pearance of hybrid synchronization. 

The net effect is that for a strong enough coupling 
both signals oscillate with the same frequency, which 
is different from the free frequencies. As found in the 
case of equal initial preparations, the fact that a remains 
mostly condensed makes a fully classical description of 
the complete system plausible. Indeed, the two frequen¬ 
cies remaining in the coupled case are reproduced by the 
classical equations [43] . In summary, for different prepa¬ 
ration of initial imbalance, it is found that the frequency 
locking with a single prominent peak appearing as shown 
in Fig. is replaced by frequency locking of more com¬ 
plex dynamics featuring several spectral components dis¬ 
played in Fig. El This result reminds us of what happens 
for two coupled classical systems, in which the coupling 
will induce measure synchronization (MS) jlHI • In classi¬ 
cal MS, the coupled dynamics will exhibit quasi-periodic 
motions, such that the Fourier analysis of Za and zi, shows 
many peaks rather than one. And in this very special 
case, the hybrid synchronization is accompanied by MS 
in the combined dynamics. We emphasize that in general 
these two phenomena do not need to arise together. 

The overall picture as Aab is varied, see Fig. is 


similar to the case of equal initial average population 
imbalance (see Fig. [^. The effect is slightly degraded, 
finding a lower condensed fraction for similar values of 
Aab in tbe case of different initial average population im¬ 
balances. We have also considered different choices of 
the nonlinearity, i.e., A^ = 2, Af, = 1. This has a similar 
effect as the different preparation of initial imbalances. 

Up to now we have only considered hybrid synchro¬ 
nization around stable phase-space points with </> = 0. In 
this case, the classical description of the Josephson junc¬ 
tion |26] shows a single stable minimum for z = 0 with 
repulsive interactions. A similar single solution, non- 
bifurcated, is found for initial preparations (pa = (j>b = 
if (Aa < l,Ah < 1). In this case, we find similar re¬ 
sults as those reported above. A more involved situation 
is found if we consider a bifurcated region of the phase 
space of each individual Josephson junction, for instance, 
(p = TT and A > 1. In this case, the classical description 
of the junction predicts a self-trapped regime |27|. To 
illustrate this dynamical regime, we have considered the 
initial condition {za){0) = Zb{0) = 0.4 and (pa = <Pb = tJ", 
with Aa = A;, = 1.2, such that the classical description 
of each junction (uncoupled) would predict a self-trapped 
regime. In the non-coupled case, see Fig. 12 (a), the pop¬ 
ulation imbalance of the condensed subsystem, Z{,, is fully 
periodic and self-trapped. The population imbalance of 
subsystem a, (Za), features a much more complicated dy¬ 
namics with no self-trapping. Coupling both subsystems 
with Aab = 0.2, Fig. 12 (b), both signals are found to 
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FIG. 11: (Color online) Absolute value of the frequency spec¬ 
trum z{u}) (Fourier transform of the average population im¬ 
balances), for different values of the coupling A;,;,. Here we 
consider different initial imbalance in the two subsystems, 
(2a)(0) = 0.4, and 2^(0) = 0.1. The frequency spectrum of the 
quantum result, ( 2 a)(a;) is given in the left panels. Zb(oj) is 
depicted in right panels. The frequency spectra are obtained 
from time series up to Tiuax = IdOtRabi. 



FIG. 12: (Color online) Evolution of {za){t) (dashed red), 
and Zb{t) (solid black) for two different values of the coupling 
(a) Aai, = 0 and (b) Aab ~ 0.2. The initial conditions are 
(2:a)(0) = 25 ( 0 ) = 0.4 and cj)a = 4>b = for the yr-phase mode, 
with Aa = Af, = 1.2, and Na = Nb = 30. 


be much more correlated. The time correlation coeffi¬ 
cient C (2 ) Zb goes from ~ 0 for the uncoupled case, see 
Fig. [TaKb), to a value close to 0.8 for the coupled case 
(Aab = 0.2). Simultaneously, the degree of condensa¬ 
tion, similarly to the case of the 0-phase mode discussed 
above, is found to increase with the coupling between the 
two systems, although the increase is less notable than in 
the zero phase case. Furthermore, we have tried different 
initial conditions for the 7r-phase mode in the bifurcated 
region of the classical phase space, and find out that due 
to the instability associated with the bifurcation [35] , the 
parameter space is much smaller compared with 0-phase 




FIG. 13: (Golor online) Evolution of the condensed frac¬ 
tion of the a system, ni (a), and time correlation coefficient 
C(za,),zb (b), for the values considered in Fig. 
condensed fraction, ni(0) = 1. 
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The initial 


mode in achieving hybrid synchronization. 


3. Initially squeezed and Fock states 

In all previous results the initial state considered was 
a condensed many-body state, i.e. condensed fraction 
ni = 1. Thus, the effect we have described up to now 
is how by coupling the quantum to a condensed system 
the condensed fraction of the quantum state was found 
to get closer to 1. For that case, the coupling to the 
condensed state was thus helping the quantum system to 
remain coherent during the time evolution. 

In this section we broaden the set of initial states to 
consider squeezed and Fock states. Squeezed states are 
particularly useful as they can be used to improve the 
efficiency of interferometers made with ultracold atomic 
systems UHlIIlllI]. In brief what we find is that cou¬ 
pling the quantum system to the condensed one has a 
similar effect as what was described before, i.e. the cou¬ 
pling mostly prevents the dephasing and thus makes the 
condensed fraction of the quantum system remain ap¬ 
proximately constant with time. 

In Fig. we consider similar conditions as in Fig. 
but with squeezed and Fock initial states. The squeezed 

_ (fc — - -— 

initial state is built as, Ck = e 2^2 where 

sets the value of the population imbalance and a sets 
the squeezing of t he sta te. The coherent state considered 
above has a = \JNa/2. A smaller value of a provides a 
squeezed initial state. In the figure we have taken tr = 1. 
The picture is very similar to the case of an initial coher- 
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FIG. 14: (Color online) Evolution of the condensed fraction of 
the subsystem as a function of time for three different values of 
hab- In the upper panel the initial state is a squeezed state, 
with an initial condensed fraction ni(0) = 0.9216. In the 
lower panel, the initial state is the Fock state, |21,9), which 
has an average population imbalance of 0.4 and an initial 
condensed fraction of 21/30 = 0.7. All other parameters are 
the same as in Fig. 


ent state. Increasing the coupling between the condensed 
and quantum subsystems the condensed fraction is seen 
to remain closer to its initial value (not 1 in this case). 

For an initial Fock state, the behavior is similar and the 
condensed fraction remains closer to its initial value for 
large enough couplings (Fig.[^ lower panel). The initial 
state is in this case the Fock state |21,9)(cfc = 
with (za)(0) = 0.4. 

Finally, in Fig. 15 we present the time correlation for 
both the squeezed and Fock states used in Fig. M For 
the squeezed case, the picture is similar to the case of 
an initial coherent preparation. The two subsystems get 
correlated for Kab > 0.01. In the Fock regime, the picture 
clearly degrades, and although a certain synchronization 
is found, it does not abide in time. 


IV. CONCLUSION 

We have considered the coupled dynamics of two ultra¬ 
cold atomic clouds, one of which is assumed to be Bose- 
Einstein condensed during the evolution. Our main find¬ 
ing is that by increasing the coupling between the two 
subsystems two net effects take place: (1) the dephas¬ 
ing produced by atom-atom interactions in the non con¬ 
densed subsystem is found to decrease as the coupling is 
increased, and (2) the coherent oscillations of both sub- 


FIG. 15: (Color online) Time correlation coefficient 
defined in Eq. ([^ versus the coupling strength Kab for the 
squeezed (upper panel) and Fock (lower panel) cases consid¬ 
ered in Fig. 


systems are found to synchronize. This phase-locking 
phenomenon is characterized by studying the evolution 
of the average population imbalance of each subsystem 
under different conditions. When synchronization ap¬ 
pears, the a state is prevented to fragment and remains 
Bose-Einstein condensed: this allows a comparison of the 
reported phase-locking with MS within a classical de¬ 
scription. Even if the role of the condensate is dominant, 
preventing the ultracold cloud from losing coherence, the 
reported synchronization differs from entrainment. The 
ultracold cloud is indeed driven by the condensate but 
the latter is also influenced by the feedback of the cloud 
and both systems evolve towards a different oscillatory 
dynamics. Synchronization is therefore hybrid (between 
a condensate driving the ultracold cloud to remain coher¬ 
ent) but mutual, being the dynamics of both clouds de¬ 
termined by the reciprocal coupling and in spite of their 
different regime. 

Our results are of relevance for future applications of 
bimodal quantum many-body systems. In particular, 
since the dephasing arising from the atom-atom inter¬ 
actions is found to disappear for large enough coupling, 
we have a way to prevent quantum many-body systems 
from dephasing. Thus, the relevant properties stored in 
the system, such as a large squeezing parameter or a 
large degree of condensation, are preserved during the 
time evolution if the system is coupled to a condensed 
one. The hybrid synchronization described, which ap¬ 
pears together with the coherent evolution, can be used 
as an observable control parameter for the phase coherent 
evolution. 




















11 


Acknowledgments 

This work was supported by China Scholarship Coun¬ 
cil, the National Natural Science Foundation of China 
(No. 11104217, No. 11205121, and No. 11402199). We 
acknowledge also partial financial support from the DGI 


(Spain) Grant No.FIS2011-24154, FIS2014-54672-P, and 
FIS2014-60343-P, the Generalitat de Catalunya Grant 
No. 2014SGR-401, and EU project QuProCS (Grant 
Agreement 641277). B. J-D. is supported by the Ramon 
y Cajal program. 


[1] Y. Kuramoto, Chemical Oscillations, Waves and Turbu¬ 
lence (Springer, Berlin) (1984). 

[2] Z. Neda, E. Ravasz, Y. Brechet, T. Vicsek, and A.-L. 
Barabasi, Nature 403, 849 (2000); D. J. Watts and S.H. 
Strogatz, Nature 393 440 (1998). 

[3] A. Arenas, A. Diaz-Guilera, J. Kurths, Y. Moreno, and 
C. Zhou, Phys. Rep. 469 (3), 93 (2008). 

[4] A. Pikovsky, H. Rosenblum, and J. Kurths, Synchroniza¬ 
tion. A Universal Concept in Nonlinear Sciences (Cam¬ 
bridge University Press, Cambridge, England) (2001). 

[5] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 64, 
821 (1990). S. C. Manrubia, A. S. Mikhailov, and D. H. 
Zanette, Emergence of Dynamical Order. Synchroniza¬ 
tion Phenomena in Complex Systems Lecture Notes in 
Complex Systems (World Scientic Publishing Co., Sin¬ 
gapore) (2004). 

[6] G. L. Giorgi, F. Galve, G. Manzano, P. Colet, and R. 
Zambrini, Phys. Rev. A 85, 052101 (2012). 

[7] G. Heinrich, M. Ludwig, J. Qian, B. Kubala, and F. 
Marquardt, Phys. Rev. Lett. 107, 043603 (2011); M. J. 
Moeckel, D. R. Southworth, E. M. Weig, and F. Mar¬ 
quardt, New J. Phys. 16, 043009 (2014). 

[8] 1. Goychuk, J. Casado-Pascual, M. Morillo, J. Lehmann, 
and P. Hanggi, Phys. Rev. Lett. 97, 210601 (2006); O. 
V. Zhirov and D. L. Shepelyansky, Phys. Rev. Lett. 100, 
014101 (2008); O. Astafiev, K. Inomata, A. O. Niskanen, 
T. Yamamoto, Yu. A. Pashkin, Y. Nakamura, and J. S. 
Tsai, Nature (London) 449, 588 (2007). 

[9] P. P. Orth, D. Roosen, W. Hofstetter, and K. LeHur, 
Phys. Rev. B 82, 144423 (2010); G. L. Giorgi, F. Plas- 
tina, G. Francica, and R. Zambrini, Phys. Rev. A 88, 
042115 (2013); M. R. Hush, W. B. Li, S. Genway, 
1. Lesanovsky, and A. D. Armour, Phys. Rev. A 91, 
061401(R) (2015). 

[10] T. E. Lee and M. C. Cross, Phys. Rev. A 88, 013834 
(2013). 

[11] G. Manzano, F. Galve, G. L. Giorgi, E. Hernandez- 
Garcia, and R. Zambrini, Sci. Rep. 3, 1439 (2013). 

[12] A. Mari, A. Farace, N. Didier, V. Giovannetti, and R. 
Fazio, Phys. Rev. Lett. Ill, 103605 (2013); V. Ameri, 
M. Eghbali-Arani, A. Mari, A. Farace, F. Kheirandish, 
V. Giovannetti, and R. Fazio, Phys. Rev. A 91, 012301 
(2015). 

[13] M. Lewenstein, A. Sanpera, and V. Ahufinger, Ultracold 
Atoms in Optical Lattices: Simulating quantum many- 
body systems, Oxford University Press (2013). 

[14] A. D. Cronin, J. Schmiedmayer, and D. E. Pritchard, 
Rev. Mod. Phys. 81, 1051 (2009). 

[15] J.-F. Schaff, T. Langen, and J. Schmiedmayer, Riv. 
Nuovo Cimmento 37, 509 (2014). 

[16] V. Giovnnetti, S. Lloyd, and L. Maccone, Science 306, 
1330 (2004). 

[17] J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and M. K. 


Oberthaler, Nature 455, 1216 (2008). 

[18] C. Gross, T. Zibold, E. Nicklas, J. Esteve, and M. K. 
Oberthaler, Nature 464, 1165 (2010). 

[19] M. F. Riedel, P. Bohi, Y. Li, T. W. Hansch, A. Sinatra, 
and P. Treutlein, Nature 464, 1170 (2010). 

[20] J. G. Bohnet, K. C. Cox, M. A. Norcia, J. M. Weiner, 
Z. Chen, and J. K. Thompson, Nat. Photonics 8, 731 
(2014). 

[21] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

[22] B. Julia-Diaz, T. Zibold, M. K. Oberthaler, M. Mele- 
Messeguer, J. Martorell, and A. Polls, Phys. Rev. A 86, 
023615 (2012). 

[23] P. Bouyer and M. A. Kasevich, Phys. Rev. A 56, R1083 
(1997); C. Orzel, A. K. Tuchman, M. L. Fenselau, M. 
Yasuda, and M. A. Kasevich, Science 291, 2386 (2001); 
W. Li, A. K. Tuchman, H. C. Chien, and M. A. Kasevich, 
Phys. Rev. Lett. 98, 040402 (2007). 

[24] M. Lewenstein and L. You, Phys. Rev. Lett. 77, 3489 
(1996); E. M. Wright, D. F. Walls, and J. C. Garrison, 
ibid 77, 2158 (1996); A. Imamoglu, M. Lewenstein and L. 
You, ibid 78, 2511 (1997); J. Javanainen and M. Wilkens, 
ibid 78, 4675 (1997); Y. Castin and J. Dalibard, Phys. 
Rev. A 55, 4330 (1997); C. K. Law, H. Pu, N. P. Bigelow 
and J. H. Eberly, ibid 58, 531 (1998). 

[25] T. Berrada, S. van Frank, R. Biicker, T. Schumm, J.-F. 
Schaff, J Schmiedmayer, Nat. Commun. 4, 2077 (2013). 

[26] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 
Phys. Rev. Lett. 79, 4950 (1997). 

[27] S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, 
Phys. Rev. A 59, 620 (1999). 

[28] B. Julia-Diaz, D. Dagnino, M. Lewenstein, J. Martorell, 
and A. Polls, Phys. Rev. A 81, 023615 (2010). 

[29] E. Padmanaban, S. Boccaletti, and S. K. Dana, Phys. 
Rev. E 91, 022920 (2015). 

[30] X. F. Li, W. Pan, B. Luo, and D. Ma, J. Lightwave Tech. 
24, 4936 (2006). 

[31] R. Gati and M. K. Oberthaler, J. Phys. B.: At. Mol. Opt. 
Phys. 40, R61 (2007). 

[32] T. Zibold, E. Nicklas, C. Gross, and M. K. Oberthaler, 
Phys. Rev. Lett. 105, 204101 (2010). 

[33] M. Abbarchi, A. Amo, V. G. Sala, D. D. Solnyshkov, H. 
Flayac, L. Ferrier, 1. Sagnes, E. Galopin, A. Lemaitre, G. 
Malpuech, and J. Bloch, Nature Physics 9, 275 (2013). 

[34] S. Kumar and C. L. Mehta, Phys. Rev. A 21, 1573 (1980); 
S. M. Barnett and P. L. Knight, J. Opt. Soc. Am. B 2, 
467 (1985). 

[35] R. Zambrini, S. M. Barnett, P. Colet, and M. San Miguel, 
Phys. Rev. A 65, 023813 (2002). 

[36] M. Heimsoth, D. Hochstuhl, C. E. Creffield, L. D. Carr, 
and F. Sols, New J. Phys. 15, 103006 (2013). 

[37] E. Boukobza, M. G. Moore, D. Cohen, and A. Vardi, 
Phys. Rev. Lett. 104, 240402 (2010). 

[38] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, 



12 


Phys. Rev. A 55, 4318 (1997). 

[39] M. Jaaskelainen, and P. Meystre, Phys. Rev. A 71, 
043603 (2005); Phys. Rev. A 73, 013602 (2006). 

[40] R. W. Robinett, Phys. Rep. 392, 1 (2004). 

[41] C. Khripkov and A. Vardi, Phys. Rev. A 89, 053629 
(2014). 

[42] S. Ashhab and C. Lobo, Phys. Rev. A. 66, 013609 (2002). 

[43] B. Julia-Diaz, M. Guilleumas, M. Lewenstein, A. Polls, 
and A. Sanpera, Phys. Rev. A 80, 023616 (2009). 

[44] I. I. Satija, R. Balakrishnan, P. Naudus, J. Reward, M. 
Edwards, and C. W. Clark, Phys. Rev. A 79, 033616 
(2009). 

[45] G. Mazzarella, M. Moratti, L. Salasnich, and F. Toigo, 


J. Phys. B: At. Mol. Opt. 43, 065303 (2010). 

[46] A. Naddeo and R. Citro, J. Phys. B: At. Mol. Opt. Phys. 
43, 135302 (2010). 

[47] M. Mele-Messeguer, B. Julia-Dlaz, M. Guilleumas, A. 
Polls, and A. Sanpera, New J. Phys. 13, 033012 (2011). 

[48] A. Hampton and D. H. Zanette, Phys. Rev. Lett. 83, 
2179 (1999). 

[49] J. Tian, H. B. Qiu, G. F. Wang, Y Chen, and L. B. Fu, 
Phys. Rev. E 88, 032906 (2013). 

[50] H. B. Qiu, B. Julia-Diaz, M. A. Garcia-March, and A. 
Polls, Phys. Rev. A 90, 033603 (2014). 



